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ABSTRACT 

A basin-scale acoustic tomography experiment was conducted in 
the northeast Pacific from May 1987 to September 1987. In this 
thesis, the stability of the forward model is analyzed. There 
are large non-linearities in the changes in travel time 
between ray paths for the four seasons. I constructed a model 
in which the change in warming in the upper 100 m of the ocean 
was due only to changes in surface solar irradiance. The value 
of the surface solar irradiance anomalies necessary to cause 
the tomography results for warming (Spiesberger and Metzger, 
1991) was computed. This value was larger than the actual 
value of surface solar irradiance anomaly which was computed 
using inputs measured by satellite (Chertock, 1989) . 
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Chapter I 
Introduction 

1.1 Motivation and Overview 

Measurement of basin scale changes in ocean temperature, at 
the surface and in the interior, is an important step in 
detecting global climate changes. Long range ocean acoustic 
tomography is an important tool for measuring the internal 
variability (Spiesberger and Metzger, 1991; Spiesberger et 
al., 1992). Satellites are an important tool in measuring 
basin-wide variations of temperature and radiation at the sea 
surface (Stewart, 1985) . This thesis compares changes in 
basin-wide temperature in the interior with changes at the 
surface. I compute the surface solar irradiance anomalies 
necessary to achieve the topographically derived warming rate 
anomalies in the upper 100 m of the ocean along the experi- 
ment ' s path . The computation is based on a model in which the 
change in warming is due only to changes in surface solar 
irradiance anomalies. This thesis also investigates the 
stability of the models used to trace sound through the ocean. 

This thesis provides an introduction to tomography and to 
surface solar irradiance at the oceans surface in chapter I. 
Chapter II describes the forward model, examines the stability 
of the model from season to season, and compares the model to 
the experimental data. Chapter III examines the estimate of 
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the surface solar irradiance anomalies and how it is computed 
Results and conclusions are given in chapter IV. 

1.2 Acoustic Tomography 

It is difficult to obtain time and space data coverage of 
the ocean by most traditional observational methods. Shipboard 
sensors can only measure one point of an observational area at 
a given time, and must move around the area taking many 
measurements to get good spatial data coverage. This can 
provide poor temporal resolution because the parameters being 
measured change during the ship's movements. Point measure- 
ments made by moored sensors provide good temporal resolution 
of one point over time, but providing sufficient areal 
coverage requires many such moorings. 

Ocean acoustic tomography was first proposed by Walter Munk 
and Carl Wunsch(1979) for measuring mesoscale (~100 km) proces- 
ses. It is a method of determining oceanic structure and 
variability through the inversion of acoustic travel-times to 
determine perturbations in sound-speed and currents and thus 
perturbations in temperature and density. Ocean acoustic 
tomography experiments are conducted by placing acoustic 
sources and receivers in the ocean. Each source transmits a 
pulse which is received at each receiver. The travel-time for 
each ray path is compared to the computed value for a model of 
known sound-speeds and the system is then "inverted" to 
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estimate the sound-speed anomalies. 

One of the big advantages of tomography is that relatively 
few instruments are recjuired to get good areal coverage of a 
large ocean area. If S sources and R receivers are placed 
around a volume and if P multipaths can be identified between 
them then there are S x R x P pieces of data per measurement. 
If point measurements are used then there is only one datum 
per instrument (S + R data) at one time. Because sound travels 
at about 1500 m s'^ in the ocean the time to collect the 
tomography data is relatively short and thus can give a three 
dimensional synoptic measurement of oceanic parameters. 

Ocean acoustic tomography was successfully used in a 1981 
acoustic tomography experiment (Cornuelle et al., 1982; 1985) 
to map sound-speed anomalies on a 300 km square area near 
Bermuda. A 1981 experiment (Spiesberger et al., 1983) measured 
Gulf Stream meanders. The 1987 experiment by Spiesberger and 
Metzger (1991) and the 1983-1984 experiment by Spiesberger et 
al . (1992) demonstrated the use of sound to observe tempera- 
ture fluctuations over ocean basins (to 4000km) . 

1.3 Surface Solar Irradiance 

The short wave radiation from the sun that is absorbed by 
the ocean, surface solar irradiance (Qj^) , is an important part 
of the heat balance of the atmosphere and the ocean. The other 
terms in the balance are latent and sensible heat fluxes and 
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outgoing long-wave radiation. The long-wave radiation from the 
ocean, Qg, can be estimated using the formula (Gill, 1992) , 

Qb = 0.985ors(0.39 - O.OSea^^) (1 - 0.673^) (1*D 

where a is Stephen's constant for a black body, Tg(°K) is 
temperature, e^ is the vapor pressure of water at the standard 
height (mb) , n^, is the fraction of sky covered by clouds, and 
the last factor, (1 - O.Sn^), is a correction for clouds. The 
total upward flux of heat, Q, from the ocean is the sum of the 
fluxes of the individual processes, namely (Gill, 1982) : ' 

Q = Qb + ( 1 . 2 ) 

v;here is the surface solar irradiance, for our purposes, 
computed by Chertock's algorithm ( 1989 ) . is the upward 

sensible heat flux, and L^E the upward latent heat flux, with 
E the evaporation rate, and the latent heat of vaporization 
of water given by 

= 2.5 X 10® J kg . (1.3) 

A simple model is used in this thesis to describe warming 
anomalies in the upper 100 m along the experiment's path. The 
model considers surface solar irradiance anomalies as the only 
source of warming anomalies. The other terms in the heat 
budget are not considered, including advection. I will use 
this model to quantify the magnitude of the surface solar 
irradiance anomaly necessary to cause the warming observed by 
tomography . 
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Chapter II 
Travel-Time Changes 



2.1 The Experiment 

The experiment used three acoustic sources and seven receiv- 
ers. The sources were deployed in a triangle, approximately 
1000 km on a side, north of Hawaii (Figure 2.1). Source (SI) 
was deployed in about 5500 m of water on a taunt mooring at 
667 m depth. It transmitted every 2 hours, 12 times a day, 
every fourth day. The pulse-like signal's acoustic frequency 
was centered at 250 Hz with an rms bandwidth of 63 Hz. Thus it 
had a pulse resolution of (63Hz)’’ « 16 ms (Spiesberger and 
Metzger, 1991) . The source level was 191 dB re 1 fiPa @lm. The 
time of the source transmission of the signal is known to 
about 1 ms with the use of a rubidium frequency standard. The 
time was known to the same accuracy at the receiver. The 
geographic positions of the sources were measured to an 
accuracy of about 10 m with the Global Positioning System. 
Mooring motion was tracked with the use of four bottom-mounted 
transponders and an interrogator on the mooring near the 
source (Liberatore, 1985) . The receivers are bottom mounted 
off the west coast of the United States. Their positions are 
known to an accuracy of about 120 m (Spiesberger and Metzger, 
1991) . 

The Acoustic path I studied extends over about 2000 km and 
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was insonified for 120 days from May 1987 to September 1987. 
It was necessary to raise the SNR (Signal to Noise Ratio) 
because the source could not transmit a pulse powerful enough 
to be heard above the background noise level . The SNR was 
increased about 36 dB when the received acoustic signal was 
complex demodulated, low-passed filtered, correlated with a 
replica (phase-only-filtering) , and averaged over four 
sequence periods (Spiesberger et al., 1989). The data were 
corrected for mooring motion and clock drift to within the 
limits noted above. The source used for this experiment is 
designated SI and the receiver is designated R1 on Figure 2.1. 
Source S3 and receiver R2 were used by Spiesberger and Metzger 
(1991) . 
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tion of the sources and receivers indicated. The acoustic path 
used in this paper is indicated with the heavy dashed line. 
From Spiesberger and Metzger, 1991. 
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2.2 Ray Tracing 



The speed of sound in the ocean depends on temperature, 
pressure, and salinity. Del Grosso's (1974) algorithm for 
sound-speed in sea water is, 

<^STP = ^^000 + ^Cp + ^CsTP ( 2 . 1 ) 



where , 

Cooo = 1402.392 

ACj. = 5.01109r- .550946xl0-^r2 + .2215153 XI O'^T^ 
Ac^ = 1.32952S + .128955X10-^52 

^Cp= .156059P+ .244998X10-^P2 _ . 883392X10-®p3 
tiCsTp = - .121xl0-'^TS + .635X10-2TP + . 2654X10-‘^r2p2 
I59xl0-^PP2 + .522X10-^TP^-.438X10-6T^P 
16 1X10 -®52p2 +.96BX10-^T^S + . A85X10'^TS^P 
34.0X1 0'^TSP 



where T is temperature in degrees Celsius, S is salinity in 
parts per thousand, and P is pressure in kilograms per square 
centimeter. Temperature and pressure have a much greater effect 
than salinity on sound-speed in most regions of the ocean. 

We will now develop a solution to the acoustic wave 
equation in terms of ray theory. The wave equation is (Tolstoy 
and Clay, 1966) 

V 2 p - = o (2.2) 
c2 at2 

where ^ is the Laplacian operator, p is the acoustic pressure 
and c is the sound-speed, a function of salinity, pressure, and 
temperature which may vary with spatial coordinates (i . e. , c = 
c(x,y,z)). For a harmonic source, with w being the angular 
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frequency, e''“^, the wave equation becomes the Helmholtz 
equation: 

V^p + k^p = 0 (2.3) 

with k being the wavenumber. Rewriting the wavenumber k as 

k = — = —■^=k^n (2.4) 

C Cq c ° 



where Cq is a constant reference sound-speed and 

n = -^ (2.5) 

c 

is the index of refraction. I write acoustic pressure in the 
form 

p = j^^ikoS{x.y.z) (2.6) 

where A is the wave amplitude, a function of (x,y,z), and k^S 
is the phase of the wave. Substituting equation (2.6) into 
equation (2.3), and collecting together the real and imaginary 
parts results in the equations 

- k^A[n^ - (V5)2] = 0 (2.7) 

and 



2VA • V5 + AV^S = 0 



( 2 . 8 ) 



If we make the assumption: 



k^A 



< 1 



(2.9) 
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then this implies the rate of variation of the wave phase, per 
wavelength, of the vertical component of wavelength is small 
(Tolstoy and Clay, 1966). Equation (2.10) can be shown to 
follow from equation (2.9) (Frisk, 1991) 



dc 

dz 



( 2 . 10 ) 



where A is the wavelength. This is a necessary, but not 
sufficient, condition for the ray acoustics approximation. It 
is also necessary that the reference point is not in the shadow 
zone or near its boundaries, and also not near caustics 
(Brekhovskikh and Lysanov, 1982) . Even simpler, the medium must 
vary slowly over an acoustic wavelength. This is the geometric 
approximation . 

Applying the condition (2.9) to equation (2.7), results in: 



(V5)2 = 



or, expanding the gradient operator: 



( 2 . 11 ) 




( 2 . 12 ) 



This is the eikonal equation, which forms the basis for ray 
theory. A physical picture of rays follows directly from the 
eikonal equation. Surfaces of constant phase (wave fronts) are 
given by S = constant, and the lines orthogonal to the wave 
fronts (VS) define the geometry of the rays (Brekhovskikh and 
Lysanov, 1982) . The rays represent the paths along which 
acoustic energy propagates. The amplitude of the rays is given 
by equation (2.8) , the transport equation. For our purposes, we 
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consider propagation in the vertical plane and let sound-speed 
be a function of depth only (i.e., c = c(z)) between sound- 
speed profile segments; a good approximation because of the 
vertically stratified nature of the ocean. Thus, if an acoustic 
signal has a wavelength that is much less than the length scale 
being examined then the geometric approximation is valid and 
ray theory can be used to model the acoustic field. 

In middle latitudes the sound-speed decreases from the 
surface to a minimum axis at ~1000m due to temperature decreas- 
ing, and then increases to the bottom due to pressure increas- 
ing in the nearly isothermal depths of the ocean (Fig 2.2). 




Figure 2.2 A typical sound-speed profile. 
Sound-speed is computed with Del Grosso ' s 
alg. using Levitus data (Levitus, 1982) . 
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Sound rays bend towards minimum sound speed according to 



Snell's law of refraction 



COS0, COS0, 

= = Constant 



( 2 . 13 ) 



where c,- is sound-speed in the layer and 0, is the angle the ray 
makes with the horizontal. 




This shows that sound rays will constantly refract toward the 
area of lower sound-speed. Because there is a minimum in sound- 
speed, at about 1 km depth at mid-equatorial latitudes, sound 
energy is bent away from the high loss regions of the surface 
and bottom. Each ray path samples a different vertical section 
of the ocean depending on the launch (departure) angle from the 
source. Travel times for ray paths are different because sound 
speeds are different along their paths. These times are given 
by 
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(2.14) 




ds 

c(s) 



where T is the travel time along the ray path for a particular 
ray, p is the ray path, ds is the arc length, and c(s) is the 
sound-speed at each point along the path. The ray paths are 
such that only a few rays actually connect the source and the 
receiver; these are defined as eigen-rays. 



2.3 The Ray Tracing Progrcua 

In this project, we used the range dependent ray tracing 
program MPP (Multiple Profile ray tracing Program) developed by 
C. W. Spofford. The sound-speed field is linearly interpolated 
in both depth and range in specified triangular sectors. The 
sound-speed field is continuous everywhere, but is's gradient 
is discontinuous at the triangular boundaries. The bottom 
bathymetry is represented by piece-wise linear segments. Output 
of the program includes eigen-ray arrival times and transmis- 
sion loss (calculated from geometrical spreading and losses due 
to boundary reflections) , along with a history of the eigen-ray 
trajectories. 

The geodesic path between the source and receiver was 
computed using the WGS84 reference ellipsoid of the Earth 
(Defense Mapping Agency, 1987) . We compute sound speed profiles 
from the Levitus seasonal data base (Levitus, 1982) using Del 
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Grosso's algorithm (ecjuation. 2.1) . Bathymetry points along the 
geodesic are from a digital data base (National Geophysical 
Data Center, 1987) . Using a flat-Earth coordinate transforma- 
tion (Ben-Menahem and Sarva, 1981) , rays are traced in a 
Cartesian coordinate system. 

MPP computes two types of rays (Spiesberger et al., 1991). 
The first type, the geometric ray, passes through both the 
source and receiver. The second type, the diffracted arrival 
(Brown, 1982; Pierce, 1981), does not pass through the receiv- 
er, but through points near the receiver. Because propagation 
is at finite frequency, energy leaks into the receiver by 
diffraction as an exponentially decreasing wave. 
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2.4 Linearizing The Forward Problem 



The travel-time of the ray along path is given by 



T, = / 

Jrj 



ds 



c(s) + U o T 



(2.15) 



where ds is a differential arc-length along the path T-, c(s) 
is the sound-speed, u is the current component along the ray, 
° is a dot product, and r is a unit vector tangent to the ray. 
The travel time of a given ray is dependent upon the path 
length, sound-speed, and current velocity along the ray path. 
Variations in sound-speed and current will lead to deviations 
in travel-time and a change in the ray path. Hamilton et al., 
(1980) showed that there is a small change in travel-time 
associated with this change in path length if sound speed and 
current variations are small. However, the change in path 
length leads to a different sampling of the oceanic medium. 

Now we will examine the relative size of the terms in the 
denominator of the integrand of equation (2.15). A typical 
value for current speed is u = 0.1 m/s and a typical sound- 



speed is c = 1500 m/s, so 



= 0 ( 10 *'*) < 1 . 



Typical values for the vertical shear of current and sound 



speed are = Q • 1^/g = 0(10'^) and 

dz 100/n 



^ . -5i2Z® =o<io-^) , 

dz 100m 



so is usually one order of magnitude larger than . 

The refraction of rays is dominated by the sound-speed gradi- 
ent, and the current can be ignored in ray tracing, to first 
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order. 



The path of the ray in equation (2 . 15) is a nonlinear 
function of the sound-speed field, so we linearize about a 
reference state, Cp, as in: 

c(X) = Co(X) ^c(x) (2.16) 

where Sc is the perturbation from the reference field. Also 

|6c(x)| < c^(X) (2.17) 



Substituting (2.16) into the travel-time integral (2.15) gives 

" fr, c^iS) + 6c{S) + u • T ■ (2.18) 

Expanding the integrand, keeping only the leading order terms, 
and dropping the current velocity term yields: 



„ f dS _ r 6c(g) dS 

^"Jr.cjS) Jr, diS) 



(2.19) 



Thus the travel-time perturbation can be expressed as 



r bc{S) dS 
c|(S) 



( 2 . 20 ) 



Equation (2.20) is linear in the sound-speed perturbation, 6c, 
if Fj is assumed not to change as the sound-speed field 
perturbs. The forward problem is finding the 6Tj when 6c is 
known . 

Sound-speed is directly proportional to temperature (6) , 
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with an empirical relationship given by: (Del Grosso, 1974) 



= a60 
Co 



( 2 . 21 ) 



where a = 3.2 X 10"^. Thus equation (2.20) can be considered a 
linear relationship between the travel-time perturbations and 
perturbations in temperature. 



2.5 Identification of Eigen- rays 



To identify eigen-rays for the spring and summer seasons I 
computed a daily average of the 12 arrival records for each day 
data was available as prescribed by Spiesberger et al. , (1980) . 
I compared the output of MPP for spring and summer seasons to 
the daily averages. The output of the MPP ray tracing program 
is shown in Figures 2.4a and 2.5a. The Levitus data base for 
the spring and summer seasons is used to generate sound speed 
profiles with Del Grosso's algorithm (equation. 2.1). The 
Levitus data base is divided into the four seasons according to 
Table 2.1. 
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Table 2.1 Levitus data base seasons 



Season 


Calendar Months 


Calendar Days 


Spring 


May - June 


121 - 212 


Summer 


August - October 


213 - 304 


Winter 


November - January 


305 - 031 


Fall 


February - April 


032 - 120 



Daily averages for Days 189 and 257 are shown in Figures 
2.4b and 2.5b. Six eigen-rays (all geometric ray arrivals) were 
identified for the spring season. The last pair of arrivals, in 
Figure 2.4b, are predicted as diffracted arrivals from MPP 
output. Because there are a large number of diffracted rays 
from MPP output grouped together near the peak arrival time, 
the individual diffracted rays could not be resolved with 16 ms 
resolution. The individual rays of the groups near the last two 
arrivals have different paths, sampling different parts of the 
ocean and thus are not used here for tomographic inversions. 

For the summer season seven geometric eigen-rays were 
identified. The last arrival was predicted from MPP as a 
diffracted arrival and the individual diffracted rays could not 
be resolved. The geometric rays that have larger travel times 
than the identified eigen-rays intersect the bottom at large 
angles and are not observed at the receiver. 

The daily average for day 189 is compared with the MPP spring 
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ray trace output in Figure 2.6a. Figure 2.6b shows the daily 
average for day 257 compared to the MPP spring ray trace 
output . 



s I v4 Forward geometries & difeacted rays from mpp, spring only 




Tkne(sccoods): Solid = geometric. Dashed = difiacted 
xl0*2 Daily average for slv4 for day 189 




Figure 2.4 (a. Top) MPP output for the spring Levitus data 
base. (b. Bottom) Daily average for julian day 189. X-axis is 
time(s) , Y_axis is relative magnitude and the start time is 
about 20 minutes . 
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si v4 Forward geometries & difiracted rays from mpp. Summer only 




Time(seconds): Solid = geomctricJDaAcd = difracted 
xlO*^ Daily average for slv4 for day 257 




Figure 2.5 (a. Top) MPP output for the summer Levitus data 
base. (b. Bottom) Daily average for julian day 257. The X-axis 
is time(s ) , the Y-axis is relative magnitude, and the start 
time is about 20 minutes. 
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X 1 0 Daily average for s 1 v4 for day 257 
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Figure 2.6 (a. Top) Day 257 average overlaid with MPP summer 
output (b. Bottom) Day 189 average and MPP spring output. The 
X-axis is time(s) and the Y_axis is relative magnitude. The MPP 
output was shifted 60ms. 
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Table 2.2 shows the difference in travel time (ms) between 
the observed daily average peak arrival time for days 189 and 
257 and the predicted eigen-ray travel time for the spring and 
summer seasons respectfully. Predicted times were less than 
their corresponding observed arrival times for all rays (Table 
2.2). A time difference of 0.1 second between measured and 
predicted travel times in the upper 1000m along our ray paths 
equates to less than 0.1 °C temperature difference (Spiesberger 
et al. , 1983) . 



Table 2.2 Daily average peak arrival time (measured) - ray 
travel time (predicted); for days 189 and 257, Spring and 
Summer . 



Peak/Ray 


Spring travel time 
difference (ms) 
measured - predicted 


Summer travel time 
difference (ms) 
measured - predicted 


4 


60 


9T 


2 


7« 


97 


3 


70 


101 


4 


HO 


l»i 


5 


118 


130 


6 


133 


156 


7 


Not Applicable 


159 
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2.6 Variations in the Forward Model 



I ran the MPP ray tracing program for the four seasons using 
the Levitus climatological data base of the oceans to compute 
sound speed profiles with Del Grosso's algorithm (equation. 
2.1). We identified six eigen-rays (all geometric arrivals) for 
the spring, fall, and winter seasons and seven for the summer 
season. The first 100 km of the propagation path of the six 
eigen-rays common to all four seasons are shown in Figure 
2.7(a-f) . 



30 




km. X-axis is horizontal range (km) and the Y_axis is depth (m) . 
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Figure 2.7 (e - f) Paths of last 2 eigen rays for the 1st 100 
km. The X-axis is horizontal range (km) and the Y-axis is 
depth (m) 



To compare the nonlinearity of the forward problem for the 
four seasons we follow in a similar fashion to a derivation 
found in Spiesberger (1985) . The reference travel time between 
a fixed source and receiver along the reference ray path is 
given by: 

T, (2.22) 

where c^ = Cq(x,z) is the initial sound speed. The perturbed 
travel time along the perturbed ray path r, is given by: 
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(2.23) 



Ti = f 

^ JFi Cp + 6c 

where 5c = Sc(x,z) is the perturbation in the sound speed. The 
perturbed ray path becomes the unperturbed ray path as 5c 0. 
The change in travel time between the initial and perturbed 
state is given by; 



- To = f — — 

^ ° Jr, Cp + 6c Jr, Cp 



(2.24) 



Because |5c/Cq| » .003 for mid-ocean sound speed perturbations 
we can expand T^ as: 



= /r. 



f dS _ 


f 


1 _ 






Jr, Cg + 6c 


Cp 


-L 


Co 


' J 



(2.25) 



Thus 5T^ becomes: 



6Ti = f 

Jr, c„ 



Cp 



1 - + 
^0 



fdcY 



\ ^0 J 



f ^ 

M Cp 



(2.26) 



or 



ff — - 


f —]- 


f dS 


6c 


6c^ 

+ • • • 


Cp 


Cp J ■ 




cl 





(2.27) 



This can be written as: 



6T, = 



ff — - 


f — 1 - 


f 


1 - 


+ . . . 


Cp 


Jro Cp j . 






Co 



(2.28) 



the last term of equation (2.28) can be written as 
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(2.29) 






1 - 



6 c 



= -/ dS-^ + 6T, 



1 



The first term of on the right-hand-side of equation (2.29), 



6 T, = -f dS 
Jro 



Co^ 



(2.30) 



is a linear approximation for the exact travel time change. The 
most important part of the linearizing is that r, = Fq (i.e. the 
perturbed path is the same as the reference path) . 

Defining 



5 T 3 = f - f 

Jr, Co Jr, Co 



and writing ST, from equation (2.29) as 



6 c 



6 c 



6T. = f dS— - f dS— + f 

* Jr. ^2 Jp ^2 Jr 



dS 6 c^ 



Co" •'r. 



1 Co Co (Co + 6 c) 



(2.31) 



(2.32) 



allows the exact travel time change to be written as 

6 T 1 = 6 T 2 + 6 T 3 + 6 T 4 . (2.33) 

This can be expressed as 

6 T 1 = Ti - To = 6 T 2 + iW: (2.34) 



or 



NL = (T^ - To) - 6T2 



(2.35) 



where 
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(2.36) 



NL = 6 T 3 + 6 T 4 . 

NL contains all of the nonlinear terms of the forward problem. 
If the initial ray and the perturbed ray coincide everywhere 
then NL = 0. The percentage of nonlinearity of the forward 
problem is: 



^NL = 100 



NL '' 

6tJ • 



(2.37) 



%NL is a measure of the nonlinear effects, primarily due to 
changes in the ray path, it is a measure of the errors of the 
linear approximation iTj verses the exact travel time change 

For our experiment we used the Levitus summer sound speed 
profiles and the MPP output ray paths as our reference state 
and the other seasons as our perturbed states. and , the 
travel times along the reference and perturbed paths are 
obtained from MPP output. %NL was computed using the summer as 
reference for all six of the eigen-rays for each of the other 
seasons. The results are summarized in Table 2.3. 

Most of the rays have large nonlinearity, thus iTj is not a 
good approximation to the exact travel time changes for these 
rays. There are exceptions; for example eigen-ray # 3 for two 
seasons and eigen-ray # 4 for all seasons have small nonlinear- 
ity. 
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Table 2.3 Seasonal variations in the forward model. <5T, is 
given by Eq. 2.24, 5 T 2 by Eq. 2.30, NL by Eq.2.35 and %NL is 
given by Eq. 2.37. The reference paths are computed from the 
summer time values. 



Season 


5T,(s) 


«T,(s) 


NL(s) 


%NL 




Eigen ray # 2 


Spring 


. 125 


.170 


-.045 


-35.8 


Fall 


. itc 


. 196 


-.076 


-63.6 


Winter 


.174 


.309 


-.135 


-77.1 


Eigen ray # 2 


Spring 


.139 


.195 


-.054 


-39.7 


Fall 


.133 


.231 


-.054 


-73.4 


Winter 


.192 


.362 


-.170 


-88^7 


Eigen ray # 3 


Spring 


.174 


.193 


-.017 


-10.0 


Fall 


.171 


.203 


-.032 


-18.3 


Winter 


.237 


.351 


-.114 


-47.9 


Eigen ray # 4 


Spring 


. 188 


.174 


.014 


7.39 


Fall 


.185 


.174 


.014 


7.42 


Winter 


.255 


.301 


-.046 


-17.98 


Eigem ray # 6 


Spring 


.Its 


.110 


.061 


36.4 


Fall 


. 188 


.068 


.120 


63.4 


Winter 


.276 


.147 


.120 


46.6 


Eigen ray # 6 


Spring 


.171 


.110 


.061 


35.6 


Fall 


.190 


.064 


.120 


66.3 


Winter 


.286 


. 142 


. 144 


50.4 
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Table 2.4 shows the first upper turning depths for the six 
eigen-rays for each season. 

Table 2.4: Eigen ray upper turning depth for the four seasons. 



First upper turning point (m) 


Ray # 


1 


2 


3 


4 


5 


6 


Spring 


0 


0 


0 


0 


18.5 


23.1 


Summer 


16.6 


22.1 


30.9 


32.2 


38.9 


40.8 


Fall 


0 


0 


0 


0 


0 




Winter 


0 


0 


0 


0 


0 


l=LJ 



2.7 Observational Variations in Travel Time 

Figure 2.8a is a gray scale bit plot of the observed data 
with the darker shades corresponding to stronger amplitudes. 
Peaks (arrivals) were picked from the data using a variable 
size window scheme that determined the largest signal within 
the window. The picked peaks are shown in white in Figure 2.8b. 
The daily average multipath arrival (peaks) trends are shown in 
Figure 2.9. Variations in the travel time of the daily averages 
are larger for the later arrivals than the early arrivals. 
Multipath arrivals start arriving earlier around year day 180 
(late June) with the sharpest decrease in arrival times between 
late June and late July. The later arrivals correspond to rays 
with smaller launch angles, thus sampling less of the deep and 
surface water. 
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Figure 2.8 (a.) Gray scale bit plot of received signal data. X- 
axis is yearday. Starting 141 and ending 261. Y-axis is tiine(s) 
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TIME (SECONDS) START TIME 




X-axis is yearday and the y-axis is time(s) 
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I 




at day 140. The Y-axis is in seconds and the X-axis is 1987 
yeardays . 
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Chapter III 

Surface Solar Irradiance Variations 
3.1 Satellite Data 

The four components of the surface heat flux budget are: 1) 
the net surface short wave radiation; 2) the net surface long 
wave radiation; 3) sensible heat; and 4) latent heat. Chertock 
(1989) generated a global seven year record of net solar 
irradiance at the ocean surface for the period between November 
1978 and October 1985. Her algorithm is based on a radiative 
transfer theory. Surface solar irradiance is the net radiation 
from the sun entering the surface of the ocean. The surface 
solar irradiance was computed by subtracting the sum of solar 
energy absorbed by the atmosphere plus the solar energy 
reflected by the ocean-atmosphere system from the solar energy 
at the top of the atmosphere. 

The solar energy absorbed by the atmosphere is a modeled 
quantity computed by Chertock (1989) . Albedo is the ratio of 
the radiation reflected by a body to the amount incident upon 
it. Planetary albedo is measured by satellite for this experi- 
ment. The solar energy at the top of the atmosphere is a known 
quantity calculated according to the astronomical sun-earth 
geometry . 

The planetary albedo used as input for Chertock 's algorithm 
came from the Nimbus-7 satellite. The Earth Radiation Budget 
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(ERB) experiments Wide Field Of View (WFOV) radiometer on 
Nimbus-7 provided the measurements (NASA, 1984b) . Twenty-two 
spectral channels of radiation were measured by the ERB 
experiment instruments. The infrared radiation emitted by the 
earth and the solar radiation reflected by the earth are 
measured by the four fixed sensors of the WFOV radiometer part 
of the ERB package. The footprint (ground resolution) of the 
WFOV instrument aboard Nimbus-7 is approximately 1000km (NASA, 
1984a) . The WFOV radiometer measures the reflected solar 
radiation in the spectral range between 0.2 and (NASA, 
1984b) . 

The Nimbus-7 satellite was launched on October 24, 1978 into 
a nearly circular, sun-synchronous orbit at an altitude of 
950km. The sun-synchronous orbit provides two passes per day, 
near local noon and midnight for equatorial crossings. Measure- 
ments by a satellite in a sun-synchronous orbit are always made 
at the same times for a given location (Stewart, 1985) . This 
produces a biased distribution of data over the diurnal cycle. 
If the parameter being measured varies diurnally, as planetary 
albedo does, the bias must be accounted for in the computa- 
tions. 

The diurnal variations in the albedo are quite large, a 
factor of 1.5 to 2 for clouds and a factor of 4 for ocean 
surfaces (Hucek et al., 1987). There are also diurnal varia- 
tions in cloud cover (Minnis and Harrison, 1984) that add to 
the diurnal variability of the albedo. 
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Albert Arking of the Goddard Space Flight Center designed a 
global reflectance model to convert the measured WFOV planetary 
albedo to a daily average planetary albedo. He used the Nimbus- 
3 directional reflectance models designed by Raschke et 
al. (1973) to estimate that 90% of the solar radiation reflected 
came from cloud-land reflectance and about 10% was from ocean 
reflectance (Kyle et al., 1986). Table 3.1 shows the Nimbus-3 
directional values derived by Raschke et al., (1973) and the 
Nimbus-7 ERB WFOV composite-directional reflectance model. 



Table 3.1 Normalized directional reflectance used in Nimbus-7 
ERB algorithms to calculate planetary albedo diurnal correction 
factors (Kyle et al., 1986) 



Cosine of solar 
zenith angle 


Ocean 


Cloud/ 

land 


WFOV 

Composite* 


0.95 


1.60 


1.56 


1.00 


0.85 


1.56 


1.56 


1.00 1 


0.75 


1.60 


1.60 


1.09 


0.65 


1.17 


1.18 


1.41 


0.55 


1.43 


1.28 


1.30 


0.45 


1.78 


1.37 


1.41 


0.95 


2.24 


1.47 


1.55 


0.«5 


2.79 


1.56 


1.00 


0.15 


3.40 


1.60 


1.78 


0.05 


4.00 


1.60 


1.84 



* Composite is 90% cloud/land plus 10% ocean. 



The planetary albedo data was processed using calibration 
adjustment procedures detailed in Kyle et al. (1985) . Daily mean 
planetary albedo data is computed (Kyle et al., 1986) from the 
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instantaneous measurements of the Nimbus-7 ERB experiment. 

To use the computed daily average planetary albedo for 
studies of the global radiation budget it is necessary to 
weight the measured solar flux by the daily average insolation 
at the top of the atmosphere (Hucek et al., 1987). \gxi the 
monthly mean planetary albedo, is computed from 



u 




i=l 



Where M is the number of days in the month and k,. is the albedo 
normalization constant for the i^^ day. F is the instantaneous 
calibrated measurement of the reflected solar radiation. is 
the amount of solar radiation that would have been reflected by 
the incident object if it had a perfectly reflecting diffuse 
surface. F and F„^^ are the average of F and for all the 
measurements within the field of view for a single pass of the 
satellite. <p is the latitude of the sub-point target area, X is 
the longitude of the sub-point target area, and tj is the day 
of the year. 
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3.2 Computation of Surface Solar Irradiance 



The solar irradiance at the ocean surface was computed by 
Chertock using her algorithm based on radiative transfer theory 
(Chertock, 1989; Frouin and Chertock, 1992; Chertock et al., 
1992) . Chertock developed a layered model of the atmosphere 
based on plane-parallel theory to compute surface solar 
irradiance. The top layer is a clear sky atmosphere which is 
above, and decoupled from, an effective cloud layer. The solar 
irradiance at the surface is a result of the radiation passing 
through the clear sky and cloud layers. The solar radiation 
incident at the top of the atmosphere is reduced by scattering 
and absorption through the clear sky layer by aerosols and 
molecules. The solar radiation is further reduced as it passes 
through the cloud layer by absorption and multiple reflections. 
Chertock used this algorithm (Chertock, 1989; Frouin and 
Chertock, 1992) to convert the monthly mean planetary albedo 
data provided by the satellite to surface solar irradiance. 

The largest contributors to atmospheric absorption of 
radiation in the solar spectrum are oxygen (O 2 ) , carbon dioxide 
(CO 2 ) , ozone (O 3 ) , and water vapor (H 2 O) (Wallace and Hobbs, 
1977) . The first two, O 2 and CO 2 , are well mixed in the atmo- 
sphere and can be considered to be at constant concentrations 
for radiation study purposes. However, the concentration and 
distribution of ozone and water vapor are highly variable, both 
in time and space, thus having a variable effect on the amount 
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of reflected radiation measured by satellite (Wallace and 
Hobbs, 1977) . 

Chertock used the algorithm created by Tanre et al.,(1986) 
to compute atmospheric absorption and scattering in the clear 
sky layer. This code, referred to as "5S” (Simulation of the 
Satellite Signal in the Solar Spectrum) , relates the radiative 
properties of the planetary system to the satellite measured 
solar radiance in clear atmosphere. The ”5S" code was used to 
compute the diffuse atmospheric transmittance due to scattering 
by aerosols and molecules. The "5S" code also calculated the 
transmittance due to gaseous absorption. The code also computed 
the portion of the measured flux that is from radiation 
reflected by the atmosphere and returned to space without ever 
reaching the cloud layer. The 5S code uses the following 
inputs; 1) the total water vapor and ozone amounts, and the 
vertical distributions of the ozone and water vapor; 2) the 
geometrical conditions; 3) The ground reflectance (a function 
of wavelength) ; and 4) an aerosol model specifying maritime or 
continental background and aerosol concentration (a function of 
visibility) . 
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3.3 Anomalies of Surface Solar Irradiance 



Chertock (1989) produced a seven year (November 1978 - 
October 1985) record of net solar irradiance at the ocean's 
surface using a numerical code to implement her algorithm. I 
provided her with Nimbus-7 ERB WFOV planetary albedo data 
(obtained from NASA's Climate Data System User Support Office 
at NASA Goddard Space Flight Center) covering the months of 
interest for the Pacific experiment. She used these data to 
compute monthly mean surface solar irradiance values for April 
1987 through September 1987. The surface solar irradiance data 
are on a 9 * latitude-longitude spacial grid box over the 
world's oceans. 

Climatological monthly means were computed by averaging the 
data for the reference period (November 1978 - October 1985) 
excluding El Nino years (January 1982 - December 1983) . Errors 
in the surface solar irradiance values are estimated to range 
between 10 and 20 W/m^ (Chertock et al., 1992). Most of the 
anomalies shown in Figure 3.1(a-f) are within the error limits 
and thus are not significantly different from zero. The anoma- 
lies of surface solar irradiance were computed by subtracting 
the climatology value from the monthly value. The variation 
from the monthly means for the months from April 1987 to 
September 1987 were computed by Chertock for each grid box in 
the Pacific experiment area. The values shown on Figures 3.1 (a 
- f) represent averages for the entire nine degree grid box. 
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APRIL 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W 
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Figure 3.1a April 1987. The X-axis is longitude and the Y- 

axis is latitude. 



MAY 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W m’^) 
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Figure 3.1b May 1987, Acoustic path between SI and R1 is 
indicated with a dashed line 
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JUNE 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W m*^) 
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Figure 3.1c June 1987, The acoustic path analyzed by Spies- 
berger and Metzger (1991) S3 - R2 is indicated with a solid 

line. 



JULY 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W m’^) 
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Figure 3. Id July 1987 
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AUGUST 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W 
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Figure 3.1e August 1987 



SEPTEMBER 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W 
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Figure 3. If September 1987 
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To examine the possible effects of changes in surface solar 
irradiance on acoustic travel times I will look at heating 
rates in the upper 100 m as Spiesberger and Metzger (1991) did. 
I will use a model in which the only parameter affecting the 
heating rate in the upper 100 m is changes in surface solar 
irradiance. The model does not consider changes in other heat 
budget terms (section 1.3) nor does it consider advection. 

To look at changes in warming, I first determined the 
average surface solar irradiance anomaly along the acoustic 
path during each month of the experiment (the dashed line in 
Figure 3.1(b)). I also determined it for the acoustic path 
analyzed by Spiesberger and Metzger (1991) (the solid line in 
Figure 2.1 and in Figure 3.1(c)). The averages were determined 
by weighing the value in each grid box by the percentage of the 
path that went through the grid box. The average monthly 
surface solar irradiance anomalies along both paths are listed 
in Table 3.2. Note that the magnitude of the anomalies are 
primarily within the error limits of the surface solar irradi- 
ance data noted previously (10 W/m^ to 20 W/m^) . 



Table 3.2 The average surface solar irradiance anomalies (w/ 
m^) along the acoustic paths for the sections Sl-Rl and S3-R2. 





May 


June 


July 


August 


Sept. 


SI - R1 


-14.5 


-17.8 


-19.4 


-24.4 


-17.2 


S3 - R2 


-14.2 


-21.7 


-21.3 


-24.3 


-17.7 



I computed the change in average monthly surface solar 
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irradiance anomaly from one month to the next, <SA. For example, 
for August to September for S3 - R2 path, 5A = -17.7 - (- 

24.3)W/m^ = 6.6 W/m^ (Table 3.2). The values of <SA are between 



1.6 and 7.2 W/m^. The standard deviation of <SA is y]2 * standard 

deviation of each monthly anomaly; namely ,/2 *20 W/m^ ~ 29 . Thus 
values of the differences of the surface solar irradiance 

anomalies are not significantly different than zero. Based on 

this model, there is no significant change in the warming rate 

from the historical mean if surface solar irradiance anomalies 

is the only cause of changes in warming. 

Spiesberger and Metzger (1991) used tomography to show that 
the warming between August 14 and September 14, 1987 was less 
than normal by 0. 2 “C±0. 056°C in the upper 100 m along the path 
from S3 to R2 (Figure 2.1). If I assume that this change in 
warming was due only to changes in surface solar irradiance, I 
can compute the surface solar irradiance anomaly decrease over 
the month necessary to achieve this decrease in warming. 
Although not realistic, this assumption allows me to determine 
the magnitude of the surface solar irradiance anomalies 
necessary to change the warming by the above amount. Using the 
formula from Wyrtki and Uhrich (1982) to compute 5T, the change 
in warming (°C) is: 



6r 



PC^D 



( 3 . 2 ) 



where p is the density of water (10^ gm/m^) , Cp is the specific 



52 



heat capacity at constant pressure (4.2 J/gm ‘C) , D is depth 
(100 m) , and 6H is the change in heating (J/m^) . For 6T = 
-0.2°C ± 0.056°C this would result in 5H = - 8.4 x 10^ ± 

2.35 X 10^ W/m^. To get this change in 5H, I computed the 
necessary change in surface solar irradiance anomaly (5A) for 
30 days with, 

hH{J/m^) . Iday ^ 1 month (3.3) 

3600 s/hour 2Ahoui 30 day 

Equation 3.3 yields -32.4 W/m^ ± 9.1 W/m^ for a month. To get 
the results of Spiesberger and Metzger using only changes in 
the monthly surface solar irradiance would require a difference 
of -32.4 W/m^ in the monthly anomalies from one month to the 
next. 

For the above results I will consider two standard devia- 
tions away from each mean. This gives Sh from tomography of 
-32.4 W/m^ ± 18.2 W/m^ and from the surface solar irradiance 
model SA is 6.6 W/m^ ± 58 W/m^. These values are less than two 
standard deviations apart. 

I wish to restate that acoustic tomography measures a 
vertical slice of the ocean integrated along the acoustic path. 
Surface solar irradiance anomalies are a monthly average over 
a 9° grid of the ocean surface. The model considers only 
changes in surface solar irradiance and neglects changes in the 
advection and other heat budget terms. Changes in the neglected 
terms may dominate changes in warming. 
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Chapter IV 

Discussion and conclusions 

4.1 The Forward Model 

We used the Levitus climatological data base with the Del 
Grosso sound speed algorithm to generate sound speed profiles 
as input to our forward model. The MPP ray trace program was 
used as the forward model. MPP produces both geometric and 
diffracted ray arrivals as output. We identified six eigen-rays 
for the spring, fall, and winter seasons and seven for the 
summer season. 

We compared the MPP output for the spring and summer to the 
1987 experiment data. The MPP forward model did a good job of 
predicting arrival times. The sound rays take about twenty 
minutes to travel between the source and receiver. MPP predict- 
ed their arrival time within 0.1 seconds (Table 2.2). Tempera- 
ture fluctuations can account for a 0.1 s offset. 

4.2 Vcuriatioiis in the Forward Model 

The seasonal stability in the ray paths of the forward model 
was examined by computing the non-linearity in the travel time 
changes as a function of season. It was shown (Table 2.3) that 
there is a large percentage of non-linearity in the changes of 
travel time (>30%) between the seasons for most of the eigen- 
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rays using summer ray paths as the reference; the summer eigen- 
ray paths are not a good approximation for the other seasons 
for tomographic inversions. 



4.3 The Observed Acoustic Data 

The travel time decreases between year-day 140 and year-day 
260. The change in travel time over the 120 day period is 
largest for the later arrivals (Figure 2.9). The later arrivals 
show decreases in total travel time which indicate faster sound 
speeds along the ray's trajectory. The later arrivals corre- 
spond to rays with smaller launch angles, that sample less of 
the near surface waters. Multipaths closest to the surface show 
much smaller changes in travel time. The earlier arrivals show 
an increase in travel time from year-day 140 to year-day 195 
and then a decrease in travel time until the end of the 
experiment, year-day 260. An increase in travel time corre- 
sponds to slower sound speeds along the ray path, indicating 
cooler water. 
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4.4 Surface Solar Irradiance 

The solar irradiance anomalies at the oceans surface, 
computed by Chertock, along the experiment's sound path range 
from -4 W/m^ ± 20 W/m^ to -35 W/m^ ± 20 W/m^. The magnitude of 
the surface solar irradiance anomalies shown in Figures 3.1 (a- 
f) are typically within the error limits of the surface solar 
irradiance data and so the changes in warming due to from them 
have no statistical significance. 

I considered a simple model in which ocean temperature 
changes in the first 100 m of the water column are a function 
only of changes in the surface solar irradiance. The model does 
not consider changes in the other terms of the heat budget 
(section 1.3) nor does it consider advection. I computed the 
surface solar irradiance anomalies (using the model) necessary 
for the change in warming found by Spiesberger and Metzger. 
They derived the change in warming from tomographic inversions 
for the path they analyzed in the experiment area. They found 
below normal heating from August to September. 
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4.5 Conclusions 



Surface solar irradiance is only one of four terms in the 
heat budget, the other three being long-wave radiation, 
sensible heat, and latent heat. We must also consider advection 
of heat when looking at the ocean's heat budget. No direct 
relationship between surface solar irradiance anomalies and 
travel time changes can be inferred without doing a complete 
heat budget including advection effects. 

Acoustic tomography has been shown to be capable of measur- 
ing the thermal properties of the ocean on a basin scale, 
including inter-annual temperature variations (Spiesberger and 
Metzger, 1991; Spiesberger et al., 1992). A long term tomo- 
graphic experiment covering the entire ocean, along with a heat 
budget computation including the terms discussed in section 1.3 
and advection, would be important in measuring changes in 
warming of the ocean. 
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